Hydrodynamically synchronized states in active colloidal arrays 
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Colloidal particles moving in a fluid interact via the induced velocity field. The collective dynamic 
state for a class of actively forced colloids, driven by harmonic potentials via a rule that couples 
forces to configurations, to perform small oscillations around an average position, is shown by 
I \ experiment, simulation and theoretical arguments to be determined by the eigenmode structure of 

( . the coupling matrix. It is remarkable that the dynamical state can therefore be predicted from 

• the mean spatial configuration of the active colloids, or from an analysis of the fluctuations near 

' equilibrium. This has the surprising consequence that while 2 particles, or polygonal arrays of 4 or 

, more colloids, synchronize with the nearest neighbors in anti-phase, a system of 3 equally spaced 

O ' colloids synchronizes in-phase. In the absence of thermal fluctuations, the stable dynamical state is 

' predominantly formed by the eigenmode with longest relaxation time. 

o : 

be understood from the eigenmode structure of the Os- 
een tensor, which depends on the geometrical arrange- 
ment oscillators. Systems of a small number of elements 
are considered, in the presence of thermal noise: between 
three and five colloidal particles are arranged on equally 
spaced average positions on a circumference, and each 
one is maintained in a driven oscillatory tangential or- 
bit, with fixed amplitude but free phase and period. The 
main result is that the fundamental (longest lived) hy- 
drodynamic mode (an easily derivable quantity, in con- 
trast to finding full solutions of the system) dominates 
the collective motion in the driven steady state. As a 
consequence, the steady state can be strikingly different 
depending on the number of particles and their arrange- 
ment. The "dynamical motifs" observed here may guide 
the analysis of cilia coordination in complex biological 
systems. 

This system is realized experimentally with optical 
traps, with fast video feedback to impose the fixed ampli- 
tude driven oscillation. The only interaction between the 
elements is through the hydrodynamic flow arising from 
the colloid movements. Brownian Dynamics (BD) simu- 
lations, in which the hydrodynamic interaction is calcu- 
lated through Osccn's tensor [l^, are compared to the 
experimental data and can also be performed readily with 
a larger number of particles in the system. The assump- 
tions for this treatment [isj are a low Reynolds number, 
particles far relative to their diameter, and a steady flow, 
both of which are satisfied in the physical context. As a 
further consequence of hydrodynamic interaction, there 
are correlations in the Brownian fluctuations of different 
particles. The dynamics of systems where spheres are 
held by fixed harmonic potentials on the vertices of ar- 
bitrary planar regular polygons has been solved exactly 
within Oseen's description of hydrodynamics [l^l, giv- 
ing a basis from which to understand the coupling in the 
active scenario. 

Optical traps are used to confine colloidal beads within 
harmonic potentials, the swtem hardware is described 
in greater detail in (lol . [l5| . In this work, a varying 



' ' Hydrodynamic coupling at low Reynolds' number is 
an important feature in biological fiows, with key con- 
r'^ ^ • sequences in apparently diverse phenomena such as the 
Q [ motility of microorganisms [3 , circulation in the brain 
^ ■ and functioning of the ear In various biological tis- 
sues, a macroscopic number of units, for example cilia, 
display synchronized dynamics, leading to metachronal 
C ' waves [3| • Nearby cilia may beat in-phase or out of phase, 
, and indeed may be in a condition where is it possible 
^ ■ to readily switch between the two dynamical states 
O ' One outstanding question is what determines the char- 
es ■ acter of the dynamical steady state. Recent progress in 
I "hydrodynamic synchronization" is reviewed in [5|, and 
^sj , an overview of low Reynolds number (Re) fiows |6|. In 
^ ' an attempt to model (both experimentally and theoreti- 
^^O , cally) the physics of hydrodynamic synchronization, two 
■ main ideas have emerged. The first is to consider the 
coupling of two or more objects driven by a constant 
force over general (pre-defined) orbits Q- Within this 
idea of "rotor" cilia models, the phase of each rotor is 
free and there can be synchronization of different rotors 
under certain conditions. This has recently been studied 
very generally in Q , extending the case of circular driving 
_ ^ force orbits [3 • A different model consists of a "geometric 
k> ' switch" , and was proposed in ; here the force is discon- 
; I ' tinuous, and this is not described within the formalism 
I of It is not yet established which of these ideas is 
most appropriate to describe a biological scenario. The 
discontinuity of force should not be discounted a-priori, 
considering the fact that molecular motors undergo dis- 
crete attachment/detachment events which couple to the 
force generation. In previous work the geometric switch 
model was investigated for the simple case of two ac- 
tive elements, showing the robustness of hydrodynamic 
coupling in the presence of noise [loj . A linear chain of 
oscillators, following the geometric switch rule with gen- 
eral actuation forces, was studied recently by numerical 
methods, in the absence of noise [llj. 

In this Letter we show how the non-equilibrium dy- 
namical behavior of the "geometric switch" model can 
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FIG. 1. (Color online) Regular arrays of actively driven colloidal particles synchronize into steady collective 
dynamical states, (a) sketch of experiment. As the number of particles A*' is varied from N = 3 to N = 5 (and beyond, 
numerically), the circle radius R is increased to maintain constant the distance between bead centers d = 2i?sin(7r/A) = 8 ^m. 
(b) Images showing one snapshot of the system, where the particles are highlighted in red. Particle positions after 20 frames 
(green) and 40 frames (blue) are overlayed. Videos are available as SM, and show clearly the A'' = 3 system performing in-phase 
oscillations, the A = 4 with neighbors in anti-phase, and A = 5 with phase locking between neighbors. The optical tweezers 
system performs image analysis on every frame to determine particle positions, and thus implement the "geometric switch" 
condition described in the text, at rates between 200 and 300 frames/s. Each colloidal particle is effectively a phase oscillator, 
undergoing a motion bound in amplitude but free in period and phase. 



number of silica beads of radius a = 1.5 (Bangslabs) 
are trapped from below by a time-shared laser beam, 
focussed by a water immersion objective (Zeiss, Achro- 
plan IR 63x/0.90 W). A pair of acousto-optical deflectors 
(AOD) allows the positioning of the laser beam in the 
(x, y) plane with sub-nanometer precision; time-sharing 
is at a rate ~ 10^ Hz, corresponding to ncghgiblc dif- 
fusion of the beads in each laser cycle. The solvent in 
which the beads are suspended is a glycerol (Fisher, Anal- 
ysis Grade) in water (Ultrapure grade, ELGA) solution 
50% w/w, giving a nominal viscosity of 77 = 6mPas at 
20°C [16[. Experiments arc performed in a temperature 
controlled laboratory, T = 21°C. The trapping plane is 
positioned (20 ± 1) /im above the flat bottom of the sam- 
ple, in a sample volume that is around 100 /im thick. 

To realize the "geometric switch" condition, an ac- 
tive driving of each colloid is implemented here, similarly 
to [l3l, but for 3,4 or 5 particles, and driving the colloids 
on segments tangential to the ring on which they are po- 
sitioned on average, see Figure[Ua). A boundary is set at 
a pre-defined particle position, a distance ^ from the min- 
imum of the currently active optical trap. The trapped 
particle moves (on average) towards the trap minimum 
and, when it crosses the boundary, the current trap is 
switched off and the other trap, with its minimum a dis- 
tance A away, is activated. Therefore the amplitude of 
oscillations is A — 2^. This process is implemented via im- 
age analysis and feedback to the AOD for laser deflection 
in the experiments (and as a condition in the numerical 
simulations). Colloidal particles are always being driven, 
and out of equilibrium in the sense that they do not reach 
the minimum of the active trap. 

The optical trap potential is harmonic to a very good 
approximation, with stiffness k in the range 1.0 to 
2.6pN//im, depending on the number of beads trapped 
(each time, calibrated with precision ±0.2pN//im). The 



relaxation time tq = j/k is of the order of G.ls. The 
experiments have been performed with A = 2 /xm, ^ = 
0.31 ^m and d ~ 8 /.jm. The peri od of an isolated oscilla- 
tor is To = 2to log[(A-0/$] [10|, which under the exper- 
imental conditions is about 0.3s. With image acquisition 
through an AVT Marlin F-131B CMOS camera, operat- 
ing at shutter aperture time of 1.5 ms, and frame rate be- 
tween 200 and 300 fps (depending on the ring size, hence 
captured region of interest) there are multiple frames cap- 
tured within the relevant timescales tq, Tq. Since the con- 
figuration is analysed experimentally only at each frame, 
the corresponding time interval should be considered as a 
feedback time. Video is acquired for over 4 minutes, i.e. 
over 48000 frames. There is typically a transient last- 
ing around a few periods before the systems reach the 
steady state discussed below. In addition to the driven 
motion, colloids are affected by stochastic thermal fluc- 
tuations fluctuations and by the net flow induced by all 
other moving particles. Experiments arc performed in- 
creasing the number of beads N, maintaining constant 
the arc-distance between neighboring beads as shown in 
Fig^b). 

Figure[2I^a) shows that for iV = 3 the three beads move 
in phase with each other, whilst for A^ = 4 (Fig. ^c)) 
the nearest neighbors are in anti-phase. The behavior of 
A^ = 5 (Fig. mje)) is an apparently more complex phase 
locking. To understand the character of the steady state 
solution we explored the properties of the Oscen coupling 
tensor, which depends on the geometry of the mean par- 
ticle arrangement. Taking a step back, we recall that in 
passive systems the theoretical description of the fluctu- 
ations of many-bead ring systems held in fixed position 
optical traps was calculated recently Il4ll . extendingorig- 
inal ideas by Polin, Grier and Ouake [ITI [18| and jl^. A 
system of N beads held in static traps in two dimensions 
has 2N normal modes, and these can be calculated ana- 
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FIG. 2. (Color online) The active tangential motion of 
beads is strongly correlated in experiments. Shown 
here are the displacements x{t) relative to the mean position, 
for systems of A'^ = 3, 4, 5 in (a), (c), (e) respectively. The 
color code red, green, blue, cyan, magenta identifies different 
beads, anticlockwise. In (b), (c), (d) the instantaneous projec- 
tions Pm{t) m = 1,2, ...N on the eigenvectors of the coupling 
matrix show the important characteristic that the mode with 
longest relaxation time (red) dominates the displacement con- 
figuration. For A'^ = 5, as for larger odd-numbered systems, 
there are two degenerate modes (red, green), which are seen 
to alternate in amplitude. The color code indicates decreasing 
relaxation time: red, green, blue, cyan, magenta. 



lytically for configurations with high symmetry [Til . [T8| . 
The active driving of each bead along its tangent direc- 
tion can be thought as introducing a constraint for each 
bead to move only along its fixed tangential direction. 
This reduces the number of modes by half, so that there 
are N normal modes for a N bead system. They can 
be obtained by projecting the two dimensional system of 
equations of motion onto the tangent vectors. Simi- 
larly to previous analysis of passive geometries [13, [3 13 
the motion of the j*''-particle originates a force on the 



-particle 



that depends on the whole 



configuration, where H is the Oseen tensor [Ic 
leads to a system of equations: 
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FIG. 3. (Color online) The dynamic state converges to 
the analytical prediction in the limit of no Brownian 
noise. In numerical simulation it is straightforward to tune 
the level of noise. Here, the trajectories and mode projections 
of a system of 7 particles are shown for temperatures of T=1K 
(a,b), 273K (c,d) and lOOOK (e,f), with other physical param- 
eters close to the experiments and K=2pN/^m, ^ = 0.3 ^m. 
Colors match those in Figures [2] plus yellow and black used 
for beads (and modes) 6 and 7 respectively. 



in which the coupling force scales as a/d, and the force 
Fi acting on the z"'-particle is harmonic and tangent to 
the ring, with the "geometric switch" rule. 



F, 
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where ±-| is the coordinate of the bottom of the har- 



particle, and it can 

< Mtl)fjit2) > = 



monic well. The stochastic force fi(t) in Eq. [T] repre- 
sents the thermal noise on the i 
be assumed that < fi{t) >= 

2kBT {H~^)ij S{ti — [HI. tifi) is a versor tangent to 
the ring, at the position , with anti-clockwise direction. 
The system of eq. [1] can be linearized for small displace- 
ments. For symmetric arrangements the coupling ma- 
trix is particularlysimple, and the eigenvectors can be 
obtained readily [20[. The eigenvectors are the normal 
modes of the coupled system. 

Knowing analytically the normal modes of the system 
constrained to tangential motion, it is possible to decom- 
pose the dynamical steady state into projections onto 
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each of the modes, and look at the mode evolution in 
between switch events. Experimentally it is clear that 
for iV = 3 and — 4, there is a single mode which has 
very high amplitude (Fig[2t^b,d)). In contrast, for = 5 
two modes have high amplitude, and they alternate pe- 
riodically (Fig[2l^f)). This behavior is confirmed by BD 
simulations performed as in [loj , see figure in Supplemen- 
tary Materials. 

Let us recall that for every even value of N , the eigen- 
mode with longest relaxation time is described by the 
pairs of adjacent spheres moving tangentially in an anti- 
phase motion, as shown in the companion paper j20| . 
For odd N with iV > 3, the longest relaxation time 
corresponds to two equivalent modes, in which the mo- 
tion of neighbors is almost in anti-phase, but on cycling 
around the ring each particle is time-shifted by +AT (or 
—AT in the equivalent mode) relative to its neighbor. 
The case TV = 3 is unique, in that the analytical theory 
shows that the longest relaxation mode is one with all 
the beads moving in phase. In Fig [5] the mode with the 
highest amplitude is consistently the one with longest re- 
laxation time. Why is the steady state dynamics largely 
captured by a normal-mode analysis, and how does the 
longest lived mode determine the steady state dynamics? 
The hydrodynamic modes are a key step in constructing 
the solutions of the dynamical system, which are given 
for each stretch between consecutive switches by linear 
superposition of the eigenmodes, and can be obtained 
anal ytic ally for the deterministic (absence of noise) sys- 
tem [201 ■ The faster the mode, the more its amplitude has 
decayed before each geometric switch. In other words, it 
is the amplitude of the longest lived mode that "dom- 
inates" at the geometric switch, and thus enforces the 
overall character of the steady state solution. 

In the presence of increasing thermal noise, the bead 
trajectories and mode amplitudes can deviate from the 
solutions of the deterministic system, as shown in Fig. H 



but the solutions remain stable. For N = 7, like N = 5 
and all other odd > 3, there is a degenerate fundamen- 
tal mode [i^l, and the dynamic state shows an alternat- 
ing amplitude of the projections onto these two modes. 
The trajectories display a fixed phase relation between 
beads, and follow each other in the order 1, 3, 5, 7, 2, 4, 6 
in Fig. [Ha,c) , and in the equivalent sequence under time 
reversal in the simulation of Fig.[31[e). The nearest neigh- 
bors are almost in anti-phase, but delayed by the small 
interval Tj^/N. Next-nearest neighbors are almost in- 
phase, with a delay 2Tm /N , and so on, describing a prop- 
agating wave. Depending on the initial conditions, at low 
noise the system will fall into a state where intervals are 
either positive or negative going around the ring. This 
dynamical state has a propagating phase. At higher noise 
(i.e. higher temperature, or weaker coupling), whilst 
the system remains overall synchronized, the propagat- 
ing character is lost over long times because the system 
is able to flip between the two equivalent states. 

The simple actively driven oscillator experiments pre- 
sented here highlight the importance of geometry in de- 
termining the leading properties of the collective steady 
state. One may envision that in biological systems, 
which present complex disordered arrangements, and a 
vast number of oscillators, the behavior highlighted by 
these simple arrays could represent the local dynam- 
ics in tightly coupled sub-systems. In this perspective, 
the small-system patterns of behavior can be thought of 
"dynamical motifs" , linked simply to the geometrical ar- 
rangement of beating elements, and that can be analyzed 
to infer the properties of the individual oscillators. 
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FIG. 4. Supplementary Figure 1 Brownian Dynamics sim- 
ulations with Oseen coupling match the experimental 
results. With parameters and markers matching those of 
the experiments of Figure [2] in the Letter. These displace- 
ments obtained computationally show that the description of 
hydrodynamic coupling via Oseen's tensor is valid. 



